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摘要 : 为 验证 稿 类 驱 避 化 合 物 与 嗅觉 引诱 物 二 氧化 碳 存 在 缔 合 作用 , 并 研究 缔 合 作用 对 败 虫 驱 避 活性 的 影响 。 本 
研究 借助 计算 化 学 的 方法 获得 缔 合 体 和 缔 合 能 量 , 利用 Gaussian View 和 Gaussian 03W 软件 分 别 构建 和 优化 二 氧 
化 碟 、22 个 蒿 类 蚊虫 驱 避 化 合 物 以 及 它们 与 二 氧化 左 缔 合 后 的 三 维 分 子 结构 , 经 Ampac 8. 16 转化 后 ,获得 它们 
的 缔 合 能 量 。 借 助 定量 构 效 关系 计算 方法 研究 缔 合 作用 对 驱 避 活性 的 影响 , 利用 Codessa 2. 7. 10 计算 获得 驱 避 剂 
和 缔 合 体 的 各 类 结构 描述 符 ,， 从 包括 缔 合 体 结构 描述 符 及 特征 描述 符 在 内 的 各 类 结构 参数 中 筛选 显著 性 参数 ， 以 
茧 类 驱 避 化 合 物 对 白 纹 伊 蚊 Aedes albopictus 的 校正 驱 避 率 的 对 数值 为 活性 数据 , 建立 结构 描述 符 与 驱 避 活性 的 定 
量 构 效 关 系 ( quantitative structure-activity relationship ，QSAR ) 模型 。 结 果 获 得 了 22 个 槛 类 驱 避 化 合 物 与 二 氧化 碳 
缔 合 的 缔 合 能 量 , 计算 显示 它们 之 间 存 在 缔 合 作用 并 且 可 以 形成 缔 合 体 ; 获得 1 个 R 为 0.9643 的 4 参数 QSAR 
模型 , 这 4 个 参数 所 对 应 的 结构 描述 符 分 别 是 COM-WNSA-3 Weighted PNSA (PNSA3 + TMSA/1 000) [Zefirov’s 
PC], f-TerCO, -Min e-n attraction for a C-O bond, M-Max 1-electron reaction index for an O atom, M-Min ( >0.1) bond 
order of an H atom, 前 2 个 参数 分 别 为 缔 合 体 的 整体 结构 描述 符 和 碎片 特征 描述 符 。 计 算 化 学 结果 表明 , MKI 
避 化 合 物 与 二 氧化 矶 存在 缔 合 作用 , 该 缔 合 作用 对 驱 避 活性 的 影响 显著 。 
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Abstract: This study attempted to demonstrate the molecular association between terpenoid repellents and 
CO,, and investigated its influence on their repellency against mosquitoes. Association energy was 
calculated using computational chemistry software: three-dimensional structures of CO,, 22 terpenoid 
repellents and their complexes that are associated with CO, were built and optimized using Gaussian View 
and Gaussian 03W , respectively, and the association energy were obtained in Ampac 8. 16. Dependence 
of repellency on the molecular association was studied by quantitative structure-activity. relationship 
method. Significant activity-affecting parameters were screened from structural descriptors of 22 terpenoid 
repellents and their complexes, as well as descriptors of characteristic fragments of complexes, which 
were all calculated by Codessa 2. 7. 10. The quantitative structure-activity relationship ( QSAR) model 
was obtained to analyze the relationship between the structural descriptors and the logarithm value of the 
corrected repelling time against Aedes albopictus. Association energies between terpenoid repellents and 


CO, were thus calculated, and the results showed that the molecular association between them was strong 
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enough to form the complexes. A statistical QSAR model with four parameters and with R^ of 0. 9643 was 
built, and the most significant activity-affecting parameters were COM-WNSA-3 Weighted PNSA ( PNSA3 
* TMSA/1 000) | Zefirov' s PC], f-TerCO,-Min e-n attraction for a C-O bond, M-Max 1-electron 
reactivity index for an O atom, M-Min ( »0. 1) bond order of an H atom, respectively. The results of 


computational chemistry calculation show that the molecular association between terpenoid repellents and 


CO, is present, and the association can affect the repellency greatly. 
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T ri Jes TEERAA H EY El, 全球 约 有 3 000 
种 , 广泛 分 布 于 除 南 极 洲 外 的 陆地 上 。 白 纹 伊 蚊 
Aedes albopictus ( Skuse) 是 我 国 重 要 的 媒介 昆虫 , 在 
国外 被 称 为 亚洲 虎 蚊 ,， 是 腹部 有 日 斑 和 银 日 斑 的 畸 
色 或 深 褐色 蚊虫 , 俗称 花 斑 蚁 。 日 纹 伊 蚊 主要 在 日 
天 吸血 , 传播 登革热 和 乙 型 脑 炎 等 , 因而 是 害虫 防 
治 领 域 的 重要 对 象 之 一 。 目 前 使 用 化 学 杀 虫 剂 是 蚊 
虫 防 治 的 主要 手段 , 但 是 其 对 背 椎 动物 的 神经 毒害 
风险 及 蚊虫 的 抗 性 增加 风险 使 蚊虫 驱 避 剂 的 研究 和 
开发 显得 格外 重要 。 败 虫 驱 避 剂 是 一 类 本 吴 无 杀 虫 
活性 , 对 蚊虫 具有 了 驱 避 作 用 的 活性 物质 ( 王 宗 德 等 ， 
2008), 表现 出 民 好 的 安全 性 和 环境 友好 性 。 近 年 
来 , 国内 外 学 者 已 开展 了 大 量 蚊 忠 驱 避 剂 方 面 的 人 研 
2C GETEXS-MIZEZKARE, 2001; Debboun et al., 2006 ) 。 
HEZK Bc H PORTU EH TRAFE, Ka. HEM 
对 人 类 和 环境 友好 的 特点 , 因而 具有 良好 的 开发 应 
用 前 景 ( 王 宗 德 等 , 2004 ) 。 

在 大 力 开发 各 种 蚊虫 驱 避 剂 的 同时 ， 驱 避 机 理 
的 研究 也 备 受 关注 , 并 取得 了 较 大 的 进展 (Ditzen et 
al., 2008; Syed and Leal, 2008; Pellegrino et al., 
2011) 。 我 们 认为 , 驱 避 化 合 物 在 涂抹 到 皮肤 之 后 ， 
是 否 与 人 体 呼出 气体 或 皮肤 分 泌 物 中 的 引诱 物 相互 
缔 合 , 并 对 驱 避 活性 产生 影响 , 在 驱 避 机 理 研 究 中 
是 一 个 非常 值得 关注 的 问题 。 国 内 外 研究 者 开展 了 
驱 避 化 合 物 与 不 同 引 诱 物 组 合 的 活性 、 行为 学 和 电 
生理 研究 (Dogan et al., 1999; Bernier et al., 2005; 
HERE, 2006), 但 从 驱 避 化 合 物 与 引诱 物 相 互 
ARAF, 对 缔 合 作用 进行 定量 描述 ,并 进一步 研 
究 缔 合作 用 对 驱 避 活性 的 影响 , 尚未 见报 道 。 我 们 
背 助 定量 构 效 关系 (quantitative structure-activity 
relationship, QSAR) 人 研究 手段 开展 了 一 些 探索 性 的 
工作 ( 雇 圣 良 等 , 2012) , 初步 推断 驱 避 化 合 物 与 引 
户 物 之 间 存 在 缮 合作 用 , 并 对 驱 避 活性 产生 影响 。 
蚊虫 驱 避 剂 的 QSAR 研究 已 有 报道 ( Suryanarayana 
et al., 1991; Katritzky et al., 2006, 2008; Natarajan 
et al., 2008; Wang et al., 2008; García-Domenech et 


molecular association ; 


quantitative structure-activity 


al., 2010) , 但 基于 缔 合 作用 对 驱 避 活性 影响 的 害 
量 关 系 计算 研究 鲜 见 报道 。 前 期 研究 ( 雇 圣 良 等 ， 
2012) 证明 这 是 研究 缔 合 作用 的 有 效 手 段 。 

二 氧化 碳 是 非 肖 重要 的 蚊虫 引诱 物 , 它 不 仅 具 
有 一 般 引 诱 物 的 特性 , 而 且 具 有 激发 和 引导 蚁 忠 飞 
行 的 作用 ( Gillies，1980), 人 研究 二 氧化 碳 与 驱 避 剂 
的 缔 合 作用 颇具 价值 ， 因此 本 文通 过 计算 屿 类 驱 避 
化 合 物 与 二 氧化 碳 缔 合 后 的 能 量变 化 ,以 确定 缔 合 
作用 是 否 存 在 , 通过 定量 计算 缔 合 体 结构 描述 符 与 
驱 避 活性 之 间 的 关系 ,从 而 对 缔 合 作用 是 否 显著 影 
啊 驱 避 活 性 做 出 判断 。 


1 材料 和 方法 


1.1 化 合 物 及 其 来 源 

本 次 研究 样本 中 包括 22 个 着 类 蚊虫 化 合 物 ， 
其 中 20 个 化 合 物 来 目 参 考 文献 Wang et al., 
2008), 全 部 是 由 oJ uki B- 派 烯 合成 的 具有 六 
元 环 结构 的 衍生 物 。 增 加 的 2 个 化 合 物 8- 产 基 别 二 
氧 葛 疆 醇和 8-A ARKEN, EEX 
验 室 中 自行 合成 , 纯度 分 别 为 97% 和 88% (EFW 
等 ,2007 ) 。 
1.2 驱 避 活性 测定 

UE (repellency ratio, RR) 的 测定 利用 小 日 
鼠 完成 。 将 小 白鼠 (50 ~ 55 HS) 腹部 去 毛 , 在 脱 
毛 区 域 均匀 地 涂抹 浓度 为 0.16 mg/cm 的 测试 化 合 
物 , 随后 将 小 日 鼠 骏 露 在 100 头 已 羽化 4 - 5 d RK 
血 的 雌性 日 纹 伊 蚊 的 败 笼 中 , 持续 2 min, WRN 
虫 呆 咬 次 数 。 计 算 公 式 如 下 : 


_ TN - TNB 
”TN 


RR x 100 (1) 


TN 为 用 于 实验 的 蚊虫 总 数 ，TNB 为 叮咬 蚊虫 
的 总 数 。 对 照 实 验 是 用 驱 避 化 合 物 的 溶剂 乙醇 通过 
同样 的 方法 进行 测试 。 利 用 式 (1) 可 以 获得 对 照 实 
Uy UN 3H X (control experiment repellency ratio, 
CERR)。 利 用 式 (2) 可 以 获得 校正 驱 避 率 ( corrected 
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repellency ratio, CRR), 它 可 以 更 准确 地 评估 驱 避 
RR - CERR 


CRR 2100 - CERR * 


100 (2) 


1.3 定量 计算 

1.3.1 缔 合 体 构建 与 几何 优化 : 应 用 Gaussian View 
4. 1 构建 22 个 蛋 类 驱 避 剂 及 其 与 二 氧化 碳 的 模拟 缔 
合体 , 采用 Gaussian 03W 进行 结构 优化 , 计算 水 平 
为 HF/6-31G(d), 得 到 最 低能 量 时 的 优势 构象 。 
1.3.2 结构 描述 符 的 计算 : 利用 Ampac 8.16 软件 


对 以 上 计算 结果 进行 格式 转换 ,再 将 转换 格式 后 的 
所 有 缔 合 体 结构 导入 到 Codessa 2. 7. 10 软件 ,计算 
它们 的 结构 描述 符 及 缔 合 部 位 的 特定 雁 族 描述 符 。 
其 中 缔 合 部 位 的 特定 碎片 描述 符 是 指 由 在 缔 合 区 域 
可 能 对 缔 合 作用 产生 影响 的 原子 、 基 团 或 分 子 片 段 
所 组 成 的 特定 碎片 经 过 计算 产生 的 描述 符 , 图 1 显 
示 的 是 8- 羟基 别 二 氢 欧 缕 醇 甲酸 酯 与 二 氧化 碳 缔 
合 的 二 维和 三 维 示意 图 。 描 述 符 包括 6 种 类 型 : 结 
构 组 成 描述 符 、 拓 扑 描 述 符 、 几 何 描述 符 、 静 电 描 
述 符 、 量 子 化 学 描述 符 和 热力 学 描述 符 。 





图 1 8- 羟 基 别 二 氧 亡 缕 醇 甲酸 酯 与 二 氧化 碳 缔 合 的 二 维 ( A) 和 三 维 (B) 示 意图 
Fig. 1 Two- (A) and three-dimensional (B) plots of the complex composed of 8-hydroxylcarveol formate and CO, 


1.3.3 缔 合 能 量 的 计算 : 缔 合 的 定义 是 两 个 或 两 
个 以 上 相同 的 或 不 相同 的 化 学 本 体 之 间 通 过 可 逆 连 
接 而 形成 复杂 物质 的 过 程 。 它 是 一 种 依靠 较 弱 的 键 
力 ( 如 配 位 共 价 键 、 氢 键 ) 结 合 的 现象 , 不 引起 共 价 
键 的 改变 ,为 可 逆 过 程 , 受 介 质 极 性 和 体系 温度 的 
X Ule] 

在 Ampac 8. 16 计算 并 转换 格式 之 后 , 可 以 详 
细 了 解构 象 的 能 量 , 通过 计算 缔 合 前 后 驱 避 化 合 物 
分 子 加 上 引诱 化 合 物 分 子 的 能 量 与 缔 合 体 的 能 量 
差 , 可 以 获得 缔 合 能 量 的 大 小 。 计 算 公 式 如 下 : 

E(binding) =E(A) +E(R) -E(AR) (3) 

E( binding) 表示 缔 合 能 量 , EC A) 表示 引 诱 化 
合 物 分 子 的 能 量 , E(R) 表 示 驱 避 化 合 物 分 子 的 能 
E, E( AR) 表 示 缔 合体 的 能 量 。 
1.3.4 定量 构 效 关系 模型 的 建立 : 应 用 Codessa 
2.7. 10 软 件 通过 局 发 式 方法 筛选 出 显著 性 结构 参 
数 , 该 方法 首先 对 分 子 参 数 进行 共 线 性 控制 , 采用 
预 处 理 的 方式 进行 快速 得 选 排除 掉 一 些 参数 , 将 结 
构 描 述 符 按 参 数 模 型 的 相关 系数 降序 排列 , 依次 引 


入 剩余 结构 描述 符 中 与 活性 或 其 他 性 质 相 关系 数 最 
大 的 结构 描述 符 。 因 此 可 以 对 结构 描述 符 快速 第 
X6, 从 而 建立 构 效 关系 计算 模型 , 它 并 不 是 对 所 有 
可 能 的 参数 组 合 进行 考察 。 

1.3.5 最 佳 定量 构 效 关系 模型 的 确定 和 检验 : 参考 
Codessa 2. 7. 10 的 手册 和 相关 文献 (Suryanarayana et 
al.,1991 ) 对 定量 构 效 关系 模型 进行 各 种 判断 、 检 
验 , 最 后 确定 最 佳 定 量 构 效 关系 计算 模型 。 

1) 模型 中 结构 描述 符 个 数 的 确定 : 根据 多 元 线 
性 回归 的 要 求 , 样本 数 和 参数 个 数 的 关系 应 当 符 
合 : nz3(k«1)(n 为 样本 数 , k 为 最 终 模型 中 的 参 
数 个 数 )。 并 且 ， 当 回归 在 增加 参数 个 数 后 ,相关 
系数 增加 不 显著 (AR <0.02 -0.05), 即 出 现 转折 
点 (breaking point) 时 ,回归 终止 , 此 时 获得 的 模型 
为 最 佳 模型 。 

2) 内 部 检验 : 根据 内 部 检验 方法 (internal 
validation) , Œ 22 个 化 合 物 分 为 A、B、C3 组 , 并 
进行 相互 组 合 。 具 体 方法 是 按 22 个 化 合 物 的 编号 
顺序 , 将 第 1,4, 7，10.…… 等 分 为 A 组 , 2, 5,8 


p 等 分 为 B 组, 剩 下 的 为 C 组 。 将 新 组 成 的 A+ 
B, A+C, B +C 3 组 导入 Codessa 2.7.10, 将 其 进 
行 结构 描述 符 的 计算 和 模型 的 建立 , 得 到 相应 的 
R, F Hi fH. 按照 所 得 模型 计算 对 应 的 剩余 的 化 
合 物 的 活性 值 , 将 这 个 预测 值 导入 到 Codessa 2.7. 
10 中 进行 验证 , 得 到 相应 的 RR, F 和 s 值 , 进行 对 
比 评价 和 判断 。 


2 结果 与 分 析 


2.1 8- 羟基 别 二 和 氢 饥 缕 醇 及 其 丙 酸 酯 的 驱 避 活性 
实验 方法 及 研究 样本 中 的 20 个 化 合 物 对 日 纹 


BERE. 定量 计算 划 类 驱 避 化 合 物 与 二 氧化 碳 缔 合 对 其 蚊虫 驱 避 活性 的 影响 


1057 


伊 蚊 的 校正 驱 避 率 值 见 Wang 等 (2008 ) , 新 增 的 两 
个 化 合 物 对 日 纹 伊 蚊 的 校正 驱 避 率 值 见 表 1。 
2.2 ” 熙 类 驱 避 化 合 物 与 二 氧化 碳 缔 合 的 缔 合 能 量 


由 式 (3) 计算 得 到 的 缔 合 能 量 列 于 表 2， 表 中 


能 量 的 单位 为 kJ/mol。 


经 过 实验 室 计 算 , 氧 键 的 键 能 大 多 在 15 ~ 40 


kJ/mol 之 间 , 由 表 中 的 数据 可 以 看 出 ,， 驱 避 剂 与 二 


氧化 碳 之 间 的 缔 合 能 都 为 正 值 , 说 明 在 自然 状态 下 


这 些 缔 合 体 的 存在 在 热力 学 上 是 可 行 的 , 虽 数 值 不 


K, 但 考虑 到 计算 是 气态 条 件 , 即使 缔 合 能 仅 为 溶 
液 中 氢 键 键 能 的 1⁄3, 仍然 表明 驱 避 化 合 物 与 引诱 
物 之 间 存 在 缔 合 作用 。 


新 测试 的 躺 类 驱 避 化 合 物 及 其 对 白 纹 伊 蚊 的 驱 避 活性 


Table1  Terpenoid repellents newly tested and their repellency against Aedes albopictus 


*1 
驱 避 化 合 物 结构 式 
Repellents Formula 
8- 羟 基 别 二 氧 欧 缕 醇 


8-Hydroxylcarveol 


8- 产 基 别 二 氢 葛 缕 醇 丙 酸 酯 


8-Hydroxylcarveol propionate 


pe 


ty" 
J 


] 


人 
44 C, H; " 


CRR: 校正 驱 避 率 Corrected repellency rate. 下 同 The same for the following tables. 


2.3 最 佳 定量 构 效 关系 模型 的 确定 

按照 最 佳 定量 构 效 关系 模型 的 确定 方法 , B Z 
确定 模型 中 结构 描述 符 的 个 数 。 含 有 不 同 结构 描述 
符 个 数 的 定量 构 效 关系 模型 对 应 的 RR 值 及 R Í D 
表 3, 它们 之 间 的 线性 关系 见 图 2。 回归 计算 的 转 
折扣 出 现在 4 参数 时 , 在 4 参数 以 后 ,， 随 春 参 数 的 
增加 , 相关 系数 并 没有 显著 地 提高 ，AR 均 小 于 
0.03, 因此 选择 4 参数 模型 为 最 佳 模型 。 

获得 的 最 佳 QSAR 模型 , 其 R 为 0.9643 , F Jy 
114.80, s^ 为 0.0004。 模 型 的 :1 检验 值 、 截 距 及 4 
个 结构 描述 符 的 具体 情况 见 表 4。 根 据 显 著 性 检 
验 , 22 个 样本 , 显著 性 水 平 a =0.01 时 , 4 个 因素 
HJ FEN 4. 31 (9E z ERI TR, 2005) , 表 4 中 模 
型 的 下 值 远 远大 于 这 个 检验 值 。 和 和 AX 分 别 表示 
回归 相关 系数 和 标准 差 , 该 模型 根据 4 个 描述 符 的 
显著 水 平 按照 降序 排列 。 

利用 所 获得 的 模型 计算 获得 的 预测 值 与 实验 测 


校正 驱 避 率 校正 驱 避 率 对 数值 
CRR lgCRR 
72.7 1.86 
1.86 
试 值 的 对 比 关 系 见 图 3。 


1.1 






0.9 
b 
C 0.8 
£= 0.7 
0.6 一 4 一 相关 系数 R? 
— 交叉 检验 相关 系数 尺 。 
0.5 
0.4 


] 2 3 4 5 6 
描述 符 个 数 Number of discriptors 


不 同 参 数 个 数 的 模型 对 应 R 及 Rs 
Fig. 2 R^ and R5, of models with different number of parameters 


图 2 


2.4 模型 的 检验 
每 一 个 亚 组 由 其 他 两 个 亚 组 组 成 的 训练 集 来 预 
m, 而 每 一 个 训练 集 的 回归 方程 是 由 相同 的 描述 符 
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Table2 Energy of the association process of terpenoid repellents and CO, 


编号 驱 避 化 合 物 驱 避 化 合 物 能 量 SAARE 缔 合 体能 量 缔 合 能 
No. Repellents Repellent energy Attractant energy Complex energy Association energy 
4-(1- 甲 基 乙烯 基 )-1- 环 已 烯 -1- 乙 醇 乙 酸 酯 
( HE) "i " Bi — 248 957.49 -74 166.32 — 323 135.45 11. 64 
4-(1-Methylethylene) -1 -cyclohexenyl-1 -ethyl acetate 
4-(1- 甲 基 乙烯 基 )-1- 环 已 烯 -1- 乙 醇 两 酸 酯 
28) "i " i — 263 982.14 -74 166.32 — 338 154. 83 6.37 


4-(1-Methylethylene) -1-cyclohexenyl-1-ethyl propionate 


8-3 Ai 5l 一 p T 
3 MENAR RN -209 312. 08 -74 166. 32 -283 489. 56 11.16 
8-Hydroxylcarveol 


8- 羟 基 别 二 氧 锡 缕 醇 甲酸 酯 
4 SA i A -252 526.61 -74 166.32 - 326 699. 68 6.76 


8-Hydroxylearveol formate 


8-359 —  EE2ERIZ, BER 
5 EA USERS CL BUT -267 545.79 -74 166.32 -341 719.61 7.50 


8-Hydroxylearveol acetate 


8-3 ARAE RA 


6 -282 570. 08 -74 166.32 — 356 743.93 7.53 
8-Hydroxylearveol propionate 
7 BERN — 181 150.87 -74 166.32 — 255 326.35 9.16 
Menthol 
eom 5 
8 KNH -239 386. 16 -74 166.32 -313 559.87 7.39 
Menthyl acetate 
A k 
9 PARR -254 410.48 -74 166. 32 -328 584. 37 7.58 
Menthyl propionate 
Wi F Bz 
10 ubl — 190 606. 22 -74 166.32 — 264 782.37 9.84 
Nopol 
诺 基 D) 
11 十 下 于 基本 — 205 589. 60 -74 166.32 -279 7162.14 6.22 
Nopyl methyl ether 
Wi F ERR 
12 b i — 220 612.74 -74 166.32 — 294 780. 89 1.83 
Nopyl ethyl ether 
诺 基 D) 
13 i PIENE -235 641.16 -74 166. 32 -309 809. 25 1.77 
Nopyl propyl ether 
j N Bib 
14 TERU: BR — 233 825.33 -74 166.32 — 307 998. 39 6.74 


Nopyl formate 
Zi MI 


15 — 248 848.29 -74 166.32 — 323 020. 48 5.87 
Nopyl acetate 
NN 

16 patus PH — 263 872.95 -74 166.32 — 338 044. 95 5. 68 


Nopyl propionate 
2 , 3 - A S DG 


17 — 175 458. 65 -74 166.32 -249 629.11 4.14 
2 ,3-Epoxypinane 
m) er à Z 酸 本 
18 HAEE N GRH -251 669.81 -74 166. 32 -325 844.15 8.02 
Endo-isocamphanylmethannol acetate 
ni grece m UQ mens 
19 MALES Ona 一 206 694. 80 -74 166.32 — 340 867. 30 6.19 


Endo-isocamphanylmethannol propionate 


型 1- 异 获 烷 基 3- 戊 醇 乙 酸 栈 
20 W 1-8HSRSCAE-3 DUST BER -311 764.15 _74 166.32 —385 939.34 8.87 


Endo-1 -isocamphanyl-3-pentanol acetate 


型 1- 异 欧 烧 基 -3- 已 醇 乙 酸 酯 
21 P1 -EIA FECE —326 792. 59 _74 166.32 — 400 967.76 8.85 


Endo-1 -isocamphanyl-3-hexanol acetate 


型 4- 异 艾 烷 基 -3- 乙 基 2- 丁 醇 乙 酸 栈 
22 pIE A IBCEN -326 713.14 -74 166.32 —400 947.74 8.28 
Endo-1 -isocamphanyl-3 -ethyl-2-butanol acetate 
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RI 不 同 参 数 个 数 的 模型 对 应 的 R 及 R, 2.00 
Table 3 R^ and R of models with different number 1.95 
= 1.90 
of parameters s= 
E 1.85 
y MA _ N MA a 
参数 个 数 相关 系数 交叉 检验 相关 系数 E 1.80 
e 
Number of parameters R? R B 1.75 
E 1.70 
1 0.5977 0.5132 Ë 1 65 
2 0. 8585 0. 8243 总 1.60 
Q 
214 0.8991 5b 1.55 
i 02 l 1.50 
4 0. 9643 0.9311 1.50 1.60 1.70 1.80 1.90 2.00 
5 0.9771 0.9516 lgCRR 计 算 值 Calculated lgCRR 
。 9842 O ocio 图 3 模型 计算 的 校正 驱 避 率 对 数值 与 实验 值 的 线性 关系 


Fig. 3 Linear relationship between calculated lgCRR and 
experimental leCRR 


表 4 最 优 4 参数 模型 


Table 4 Best four-parameter model 


描述 符 序 号 回归 系数 回归 系数 标准 误 -检验 值 描述 符 名 称 
Descriptor number X AX t-test value Name of descriptors 
0 — 11. 1200 0.9738 -11.4191 Intercept 


COM-WNSA-3 Weighted PNSA 


! -0. 1023 0. 0054 - 13. 0090 ( PNSA3 « TMSA/1000) [ Zefirov' s PC] 

2 0. 0284 0. 0021 13.3496 f-TerCO, -Min e-n attraction for a C-O bond 

3 - 16. 1880 2.7227 -5,9457 M-Max 1-electron reactivity index for an O atom 
4 1.5278 0. 3284 4. 6528 M-Min ( »0.1) bond order of an H atom 


COM-; 基于 缔 合 体 计 算 产 生 的 描述 符 Descriptors generated from complexes; f-TerCO;-: 基于 特定 碎片 计算 产生 的 描述 符 Descriptors generated 
from specific fragments; M-: 基于 驱 避 剂 单 体 计算 产生 的 描述 符 Descriptors generated from repellent monomers. 


构建 。 用 于 预测 的 QSAR 模型 的 有 效 性 是 用 每 个 训 。 ”系数 (RE ) 来 评价 。 计 算 结果 见 表 5。 
练 集 的 留 一 法 (leave one out, LOO) 交 义 验证 相关 


RS 4 参数 模型 检验 


Table 5 Validation of four-parameter model 


化 合 物 个 数 交叉 检验 化 合 物 个 数 交叉 检验 op 
训练 集 相关 系数 "M 标准 差 测试 集 相关 系数 "" 标准 差 
Number of 相关 系数 Number of 相关 系数 
Training sets R° ) s? Testing sets R? ) s? 
compounds R, compounds R, 
A+B 15 0. 9853 0. 9528 0. 0002 C 7 0. 8678 - 0. 0013 
A+C 15 0. 9687 0. 9337 0. 0003 B 7 0.9145 - 0. 0006 
B+C 14 0. 9546 0. 9248 0. 0006 A 8 0.9375 - 0. 0010 
TOf 0. 9695 0.9371 0. 0004 TM 0. 9066 0. 0010 
Mean Mean 


由 于 检验 组 所 含有 的 化 合 物 的 数量 化 合 物 的 不 ”获得 的 QSAR 模型 具有 很 好 的 稳定 性 和 预测 能 力 ， 
K, 获得 的 RR 的 预测 值 会 比较 离散 , 但 是 从 这 3 个 “该 模型 可 以 很 好 地 评估 这 22 个 缔 合 体 与 驱 避 活性 
检验 组 之 间 仍 然 可 以 得 到 比较 显著 的 结果 , 说 明 所 ”之 间 的 关系 。 
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3 讨论 


本 研究 使 用 计算 化 学 手段 对 22 个 萌 类 蚊虫 驱 
避 化 合 物 及 其 与 二 氧化 矶 形成 的 缔 合 体 进行 结构 优 
化 、 描 述 符 计算 、 定 量 构 效 关系 建 模 ，, 最 终 得 到 了 
缔 合 的 能 量 和 对 应 于 校正 驱 避 率 对 数值 lgCRR 的 
最 佳 定 量 构 效 关系 模型 ,计算 结果 表明 昔 类 张 避 化 
合 物 与 二 氧化 碳 的 缔 合 作用 对 驱 避 活性 具有 显著 的 
影响。 

缔 合 体能 量 小 于 驱 避 化 合 物 与 引诱 物 的 总 能 
量 ,， 从 理论 上 说 明 缔 合体 在 热力 学 上 是 稳定 的 , 缔 
合作 用 的 存在 可 进一步 增加 其 稳定 性 。 后 续 的 研究 
有 竺 全面 地 计算 各 类 驱 避 化 合 物 和 多 种 主要 的 引诱 
物 之 间 的 缔 合 作用 , 并 从 行为 学 和 触角 电 生 理 实验 
加 以 验证 。 最 佳 定量 构 效 关系 模型 共 包 售 4 个 结构 
描述 符 : COM-WNSA-3 Weighted PNSA (PNSA3 + 
TMSA/1 000) | Zefirov s PC], f-TerCO,-Min e-n 
attraction for a C-O bond, M-Max 1-electron reactivity 
index for an O atom, M-Min ( »0.1) bond order of an 
H atom。 从 最 佳 定 量 构 效 关系 模型 可 以 知道 , 在 影 
啊 知 类 蚊虫 驱 避 化 合 物 活性 的 4 个 结构 描述 符 中 ， 
第 1 个 描述 符 属于 缔 合 体 , 是 缔 合 体 正 电荷 加 权 部 
分 表面 积 。 第 2 个 描述 符 属于 缔 合 体 碎片 , 是 缔 合 
体 碎 片 内 部 碳 氧 键 电子 与 原子 核 最 小 引诱 力 。 第 3 
个 描述 符 和 第 4 个 描述 符 属 于 驱 避 剂 单 体 ,分 别 是 
驱 避 剂 单 体 的 氧 原子 最 大 1- 电子 反应 指数 和 氧 原 
子 的 最 小 键 级 ( >0.1)。 第 1 个 和 第 2 AT XS EAD 
是 缔 合 的 结构 描述 符 , 是 影响 驱 避 活性 最 为 显著 的 
两 个 描述 符 。 第 3 个 和 第 4 个 描述 符 都 是 驱 避 化 合 
物 单 体 的 结构 描述 符 , 说 明了 了 驱 避 化 合 物 自 号 对 驱 
避 活 性 的 影响 同样 是 不 可 忽视 的 。 

由 上 分 析 可 知 , 缔 合 体 及 驱 避 化 合 物 单 体 对 豫 
避 活 性 都 存在 影响 。 缔 合作 用 对 驱 避 活性 的 影响 可 
能 是 多 方面 的 ,比如 驱 避 化 合 物 与 引 记 物 缔 合 可 能 
使 引诱 物 的 引诱 能 力 损 失 , unl SEXES C REA X 
虫 嗅觉 感受 右 后 产生 了 驱 避 作用 。 缔 合作 用 究 竞 是 
如 何 有 影响 驱 避 活性 的 , 还 有 得 更 为 深入 的 研究 , 这 
方面 的 深入 人 研究 对 于 驱 避 机 理 研 究 也 极 具 意义 。 


致谢 ”本 研究 得 到 江西 省 竹子 种 质 资源 与 利用 重点 
实验 室 的 条 件 支持 ,江西 农业 大 学 理学 院 2008 级 
应 用 化 学 专业 赖 忠 才 同学 参与 了 部 分 计算 工作 , 在 
此 表示 衷心 感谢 。 
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